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A MONTE CARLO PROGRAM FOR TRANSMISSION PROBABILITY CALCULATIONS 

INCLUDING MASS MOTIONS 

SUMMARY . 

V 

A Monte Carlo computer program which can be used to calculate 
the transmission probabilities for two cylindrical tubes in series 
is described, as well as the technique of adding directed mass motion 
to the random thermal motion of the molecules. The results for a 
simple straight cylindrical duct are compared to other solutions, and 
angle-of-attack effects are examined. Complex systems for which no 
adequate solutions have previously been found are easily analysed by 
this method. The Fortran instructions are listed along with a typical 
program solution. 


I. INTRODUCTION 

While the properties of rarefied gas flow have been studied in 
detail for many years, there still exists, in many respects, a real 
lack of understanding of the physics of the problem. Gross properties 
such as the transmission probability (the probability that, if a molecule 
enters a system at point 1, it will exit the system at point 2), the 
spatial distribution at, near, or in the system, etc., have been deter- 
mined for very simple systems (i.-e., cylindrical ducts, rectangular ducts, 
infinitely wide parallel plates, etc.) when the density is sufficiently 
low that the system is in the free molecular flow regime. However, the 
condition for free molecular flow is not exact since all the parameters 
which affect the flow properties are not known. The Knudsen number is 
normally defined as the rates of the mean free path (A) to the character- 
istic dimension (L) of the system. When the Knudsen number is large 
(Kn = 10 or higher), the system is said to be in free molecular flow. 

But if one is considering a cylindrical duct as a system, it is not 
known if the characteristic length is the diameter of the duct, the 
length of the duct, or some dimension resulting from a combination of 
these values. For very large mean free paths and relatively small 
length-to-diameter ratios for the duct, the exact Knudsen number is 
not important. When the mean free path of the molecule is near the 
characteristic length (10 ^K ^ .1)? the flow regime is called transition 
flow, for which no sufficient description is presently available. 



These problems are greatly magnified when one considers complex 
systems such as orifice-restricted tubes, baffles in tubes, elbows, etc., 
and flow where there is relative motion between the system and the gas. 
There is a wealth of information in the literature about this problem 
with most of the solutions anchored on the excellent work of Clausing [ l] . 
His approach and techniques have been followed with some improvements in 
accuracy. It was not until Davis in 1959 [2] that any significantly 
different approach to the problem was developed. In a series of extremely 
useful papers [2, 3, 4, 5], Davis and his co-workers made a significant 
contribution to the field of vacuum system design. Because of the ' 
obvious success of Davis 1 approach and the belief that this approach 
could be used for solutions of many problems, the writer has developed 
programs very similar to those of Davis to look at several types of 
problems. Some results of these programs have been published [6,7]. 

These programs consider complex geometrical systems, relative mass 
motion, reflection coefficients, spatial distributions, etc. Since 
many requests for information about these programs have been received, 
it was decided to present in detail the program for a simple system. 

This paper describes the Fortran computer program for the Monte Carlo 
calculations of the transmission probability for two co-axial cylinders 
connected in series with provisions for relative mass motion directed 
along the axis of the system or at any arbitrary angle with respect to 
the axis. Only the basic program will be described in detail. Modifica- 
tions which have been incorporated to determine flux distributions, 
spatial distributions, capture coefficients, etc., will be only briefly 
mentioned since these aspects can be added to the basic program as the 
user desires. The program has not been refined. It was not the intent 
of this study to develop an optimum, time-saving computer solution. A 
skilled programmer could easily reduce the computing time. However, it 
was felt that the 8 minutes it takes for the present computer solution 
on a CDC-3200 (or 2 to 4 minutes on an IBM-7094) was reasonable. Most 
of the effort has been spent on the physics of the problem. 


II. DESCRIPTION OF THE PROGRAM 


A. ASSUMPTIONS 

The following basic assumptions of the flow parameters of the 
molecules are consistent with most approaches: 

1. The molecules enter the cylindrical duct with a direction 
proportional to the cosine of the angle between the normal to the entrance 
plane and the direction of the molecule. If the molecule strikes a wall, 
it leaves the wall in a direction proportional to the cosine of the angle 
between the normal to wall at the point of collision and the direction 
of the molecule. 
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2. The mean free path is sufficiently large such that collisions 
between molecules can be ignored. (Modifications to this basic program 
are not limited by this assumption; however, no description of these 
modifications will be made in this paper). 

B. MODEL DESCRIPTION 


The model for this study is shown in Figure 1, where we see 
the two cylindrical tubes in series with a common axis. For convenience, 
the radius of tube 1 is given the value of one. The configuration to be 
studied is then described in input parameters for the following parameters. 
(The symbols enclosed by parentheses refer to the parameters in the Fortran 
program, Figure 3): 


1 . 

The 

2. 

The 

3. 

The 

4. 

The 

5. 

The 

6. 

The 


length of tube 1, (Q). 
length of tube 2, (QL2). 

radius of the entrance orifice to tube 1, (A), 
radius of the exit orifice of tube 1, (B) . 
diameter of tube 2, (CL2) . 

radius of the exit orifice of tube 2, (DL2). 


In Figure 1, tube 2 is shown to be larger in diameter than 
tube 1. The program is not restricted to this configuration and can 
be used with tube 2 having any diameter. If QL2 is zero, tube 2 is 
not considered by the program. 


C. THE MONTE CARLO PROCEDURE 


In the Monte Carlo procedures, molecular histories are generated 
for a large number of molecular paths through the system. The properties 
of each molecule are recorded depending on the parameters of interest. 
This paper will consider only the tabulations of the number of molecules 
which exit the system at either exit plane as a function of the number of 
collisions with the walls (both the tube walls and the orifice plates). 

A flow diagram for the computer program is shown in Figure 2; 
the Fortran print-out for the program is shown in Figure 3. The history 
for a molecule begins with instruction 2010. (The program instructions 
preceding number 2010 merely performs bookkeeping functions, presets 
counters, reads the input parameters, etc.) The initial coordinates of 
the molecule on the entrance plane are given by' (with the axis of the 
cylinders being the positive X axis) 
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Yi = 1 - 2 R x 
Z x = 1 - 2R 2 , 


where and R 2 are random numbers. (There are numerous methods to 
calculate pseudo-random numbers in the literature.) These coordinates 
can fall anywhere within the square bounded by y = ±A , and z = ±A. 

The next instruction selects only those points within the circle 

Yi 2 + Zi 2 = A 2 . Instructions 2011 through 2142 calculate a cosine-bias 

angle . 

If mass motion is to be considered, the instructions from 2555 
through 2141 are used to calculate the random thermal speed of the 
molecule. This is explained more fully in the next section. Direction 
cosines for the molecular trajectory are determined from 2556 through 
2020. Instruction 2024 tes-ts to see if the molecule hits the exit 
plane of the tube. If it does not hit this plane, it must hit the 
wall of tube between X = 0 and X = L(Q) . The point of collision is 
calculated by instructions 3026 through 3035. At this point, the 
molecule is reflected diffusely from the tube wall by the instructions 
3039 up to 3045. The instruction before 3045 senses Ul to see which 
direction the molecule will go. Instructions 3045 to 3052 and 3060 to 
3063 tests to see if the molecule hits the end planes of the cylinder. 

If it does not, the x coordinate of the point of collision with the 
tube wall is calculated by 3068 to -3642. The remainder of the program 
repeats these tests for molecules in tube 2. 

D. SPEED DISTRIBUTION 

Although the procedure for combining the relative mass motion 
to the random thermal motion is straightforward, it needs some explanation. 
In the program, the position of a molecule at the entrance plane of the 
system is determined, and the direction cosines for its motion considering 
only the cosine bias are calculated. For no mass motion, the molecules 
whose direction cosine Ul, because of its thermal velocity, is negative 
will not enter the system. This is not necessarily so when there is mass 
motion. The magnitude of the thermal speed of the molecule is chosen 
by the instructions starting at 2555 in the program. Figure 5 presents 
a histogram of the thermal speed as selected by this procedure. From 
this distribution, one can determine the most probable speed to be 1.414, 
the root-mean-square speed to be 1.777, and the average or mean speed to 
be 1.596. Also shown in Figure 5 is the theoretical Maxwellian-Boltzman 
distribution. One can see that this calculation satisfactorily approximates 
the Maxwellian distribution. 

Using the direction cosines, one can then determine the magnitude of 
each component of the speed. Mass motion can then be added to the thermal 
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speed and new direction cosines for the actual path of the molecule 
relative to the moving tube. One small consideration must be made, 
however. This consideration is concerned with the relationship of 
the thermal speed of the molecules entering the system to the thermal 
speed of the molecules in the reservoir. From the kinetic theory of 
gases, one knows that the number of molecules striking the small area 
is given by 


\ v 


where V is the speed of the molecule and 


d% - 



e 3 v 2 


-P 2 v 2 


dV 


r = J2_ 
P 2kT 


The most probable speed of the molecules is found by setting the first 
derivative of function above equal to zero. 


_d_ 

dV 



P 3 V 2 


■p 2 v 2 


= 0 


From this we find that the most probable speed of the molecules leaving 
the reservoir is given by 


^exiting 
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The distribution shown in Figure 5 is for the molecules exiting the 
reservoir (or entering the tube). Since the speed ratio is defined 
as the ratio of the mass motion to the most probable motion of the 
gas molecules in the reservoir, the mass motion vector which is added 
to the thermal speed vector must be expressed in terms of the reservoir 
speed. Thus, the mass motion speed vector is given by 


1.414 S 


1.155 S, 


where S = speed ratio. 

Figure 6 shows a histogram for the speed distribution with mass 
motion. The solid curve represents the Monte Carlo solution for the 
general expression of the velocity distribution 


d 3 N v 

N 


\2jtkT/ 


2kr [ < V x-SD 2 + V y 2 + V z 2 ] 


dV x dVy dV 2 


where N = 10,000 and SI = 8 ( ^reservoir 


mp 


The root mean square speed was found to be 9.384 and V, 9.330. The 
dashed curve shows the distribution of molecules that actually enter 
the tube for the same mass motion. It is seen that, although the rms 
speed and the mean speed are approximately equal, the actual distributions 
are very different. This difference is a result of the cosine distribution 
of the entering molecules. In Figure 6, the peak on the program distri- 
bution at V = 8.5 results from the addition of the mass motion to those 
molecules at the entrance plane which were traveling opposite to the 
direction of mass motion, while the peak at V = 10.0 results from the 
addition of mass motion with the molecules which were traveling in the 
same direction as the mass motion. 


6 



Ill, RESULTS 


The transmission probabilities for cylindrical ducts at various 
speed ratios have been calculated by Pond [8], by deLeeuw and Rothe [9], 
and by Hughes [ 10] . While Pond presents his data directly as transmission 
probabilities, it is quite difficult to accurately obtain values from his 
figures. Rothe and deLeeuw present their data as a pressure ratio; some 
simple calculations are required to convert the data into transmission 
probabilities. The results of this operation are presented in Table 1 
and in Figures 7 through 9. Figures 10 through 15 present data from the 
Monte Carlo solutions with the solid line representing the values inter- 
polated from the UTIA data. It is quite noticeable that the Monte Carlo 
data for all L/A values are slightly lower than the UTIA data below speed 
ratio ot 5. For the zero speed ratio data, the Monte Carlo values are 
consistently within ±5 per cent of other solutions. For speed ratios 
up to 4, the difference between the Monte Carlo data and the UTIA data 
is between 5 and 10 per cent; however, the Monte Carlo data are always 
the smaller. No explanation of this can be given at this time. 

All the previous results are for the case where the mass motion is 
directed along the axis of the tube. In general, the mass motion will 
be at any angle with respect to the axis of the tube. Defining the 

angle of attack to be the angle between the direction of mass motion 

and the axis of the tube, Figure 16 presents the Monte Carlo data for 
tubes of L/A = 2, and speed ratios of 2, at various angles of attack. 

It is interesting to note that, at an angle of attack of 90 degrees, the 

transmission probability (0.35) is not the same as the transmission 
probability for zero speed ratio (0.51). 

In all the data presented so far, only a single tube nas been 
considered. Normally, a measuring system in free molecule flow consists 
of a gauge volume connected to the free stream by an orifice-restricted 
tube or a series of tubes. The techniques for calculating the transmission 
probability for composite systems with zero mass motion have been examined 
by several investigators [7,11,12]. To see if these same approaches could 
be used for composite systems with mass motions, Monte Carlo calculations 
were made for an orifice-restricted tube of L/A = 8, A^/A? = .5, and for 
various speed ratios. So that any differences which might appear would 
be the result of the procedure, the Monte Carlo transmission probability 
values for a tube of L/A = 4 were used in the calculations. Figure 17 
presents the results of these calculations. Deviations up to approximately 
15 per cent are observed between the calculated values and the Monte Carlo 
values. Tube lengths of L/A = 2 and 4 were also examined (although not 
presented in this report). Deviations from the theoretical values for 
these tubes were quite small, i.e., 2 to 4 per cent. The standard tech- 
niques for coupling orifices and tubes in free molecular flow may be used 
for mass motion studies. The deviation from these values increased with 
increasing tube length, but decreased with increasing mass motion. 
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For the case of two tubes with transmission probabilities Q ± and 
corresponding to tube Lx/A-l and L 2 /A 2 , respectively. Reference 7 shows 
that the total transmission probability for zero mass motion is given by 


_i_ = j_ + A i 2 _ . 

a 2 2 qi 2 a 2 2 


Assuming this equation to be valid for the case where there is mass 
motion such that Ot( L/A) becomes G(L/A,S), calculations were made for 
two values of L/A and various speed ratios and are compared with Monte 
Carlo data for the same parameters in Figure 19. There is excellent 
agreement between the calculated values and the Monte Carlo values. 

Figure 20 presents a comparison of Monte Carlo solutions for trans- 
mission probabilities for cylindrical ducts and angles of attack with 
data from Reference 10 for the same configuration. Again, there is 
excellent agreement between the two metrhods. 


IV. CONCLUSIONS 

A simple Monte Carlo computer program similar to the one developed 
by Davis has been modified to permit investigation of the influence of 
mass motion with random orientation with respect to the tube. Comparison 
of results from this program to numerical solutions of other investiga- 
tions show the agreement of the different methods to be excellent. Using 
the inherent advantages of this program, such as orifice studies, series 
of tubes, etc., the general coupling equations developed for zero mass 
flow appear to be applicable to flows with relative mass motion. 

After this comparison of numerical methods to Monte Carlo methods 
for simple tubes, a great amount of confidence may be placed in attempting 
investigations of parameters such as reflection coefficients, thermal 
and accommodation coefficients, etc. Also rather simple modifications 
to this program allow one to examine ducts of other types, i.e., 
elliptical and conical. 
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TABLE I 


TRANSMISSION PROBABILITIES FOR CYLINDRICAL DUCTS 
AT VARIOUS SPEED RATIOS 

(CALCULATED FROM UNIVERSITY OF TORONTO REPORT NO. 88. 
DECEMBER 1962* BY J. H. DELEEUW AND D. E. ROTHE) 


L/A K( S.L/A) 



S=>0 

S=0.02 

S=0. 05 

S = 0.1 

.010021 

.995015 

.995073 

.995179 

.995317 

.013338 

.993375 

.993433 

.993630 

.993760 

.017752 

.991202 

.991355 

.991456 

.991753 

.023628 

.988324 

.988477 

.988668 

.989040 

.031449 

.984520 

.984767 

.985043 

.985482 

.041859 

.979501 

.979747 

.980111 

.980788 

.055714 

.972901 

.973239 

.973775 

.974507 

.074155 

.964256 

.964684 

.965299 

.966416 

.098701 

.952988 

.953595 

.954368 

.955765 

.131371 

.938402 

.939091 

.940191 

.941927 

.174854 

.919688 

.920629 

.921947 

.924150 

.232731 

.895950 

.897040 

.898725 

.901429 

.309765 

.866283 

.867671 

.869681 

.872967 

.412298 

•829888 

.831538 

.833904 

.837829 

.548769 

.786257 

.788124 

.790854 

.795436 

.730410 

.735367 

• 737397 

.740476 

.745562 

.972176 

.677841 

.679974 

.683233 

•688666 

1.293971 

.614986 

.617159 

.620498 

.626050 

1.722267 

.548766 

.550917 

.554137 

.559610 

2.292342 

.481741 

.483769 

.486810 

.491991 

3.051106 

.416713 

.418547 

.421327 

.426000 

4.061021 

.355790 

.357391 

.359827 

.363959 

5.405230 

.299834 

.301212 

.303291 

.306845 

7.194348 

.249045 

.250213 

.251962 

.254952 

9.575701 

.203677 

.204633 

.206081 

.208559 

12.745267 

.164051 

.164837 

.166018 

.168025 

16.963960 

.130320 

.130944 

.131882 

.133499 

22.578970 

.102309 

.102799 

.103545 

* 104813 

30.052592 

.079535 

.079916 

.080503 

.081495 

40.000000 

.061367 

• 0 6 1 6 6 7 

.062114 

.062879 
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TABLE I CONTINUED 


TRANSMISSION PROBABILITIES FOR CYLINDRICAL DUCTS 
AT VARIOUS SPEED RATIOS 

(CALCULATED FROM UNIVERSITY OF TORONTO REPORT NO. 88* 
DECEMBER 1962* BY J. H. DELEEUW AND D. E. ROTHE) 


L/A 


K(S*L/A) 



S = 0. 2 

S = 0 • 3 

S = 0. 5 

S=0 . 75 

.010021 

.995636 

.995930 

.996414 

.997009 

.013338 

.994209 

.994534 

.995239 

.995994 

.017752 

.992247 

.992786 

.993668 

.994687 

.023628 

.989720 

.990392 

.991619 

.992945 

•031449 

.986405 

.987248 

.988868 

.990612 

.041859 

.981939 

.983122 

.985209 

.987524 

.055714 

.976090 

.977638 

.980400 

.983435 

.074155 

.968454 

.970379 

.974047 

.978019 

.098701 

.958367 

.960922 

.965667 

.970845 

. 131371 

.945314 

.948589 

.954683 

.961429 

.174854 

.928441 

.932566 

.940359 

.949078 

.232731 

.906855 

.912034 

.921860 

.932991 

.309765 

.879561 

.886005 

.898220 

.912209 

.412298 

.845703 

.853442 

.868328 

.885554 

.548769 

.804512 

.813570 

.831181 

.851933 

.730410 

.755763 

.765950 

.786133 

.810390 

.972176 

.699656 

.710758 

.732956 

.760292 

1.293971 

.637380 

.648912 

.672337 

.701856 

1.722267 

.570796 

.582290 

.605929 

.636333 

2.292342 

• 502600 

.513549 

.536385 

.566302 

3.051106 

.435683 

•445744 

.466882 

.494972 

4.061021 

.372523 

.381433 

.400314 

.425685 

5.405230 

.314215 

.321907 

.338272 

.360475 

7.194348 

.261150 

.267640 

.281523 

.300456 

9.575701 

.213679 

.219073 

.230612 

.246424 

12.745267 

. 172202 

. 176594 

.186016 

.198984 

16.963960 

. 136841 

.140364 

.147934 

.158386 

22.578970 

. 107458 

.110245 

.116238 

.124537 

30.052592 

.083566 

.085744 

.090446 

.096955 

40.000000 

.064482 

.066169 

.069811 

.074866 
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TABLE I CONTINUED 


TRANSMISSION PROBABILITIES FOR CYLINDRICAL DUCTS 
AT VARIOUS SPEED RATIOS 

ICALCULATED FROM UNIVERSITY OF TORONTO REPORT NO. 88* 
DECEMBER 1962» BY J. H. DELEEUW AND D. E. ROTHE) 


L/A K ( S*L/A ) 



S=1.0 

S = 1 • 5 

S = 2 • 0 

S=3 • 0 

.010021 

.997460 

.998153 

.998587 

.999066 

.013338 

.996636 

.997533 

.998118 

.998727 

.017752 

.995520 

.996726 

.997500 

.998314 

.023628 

.994043 

.995648 

.996680 

.997738 

•031449 

.992060 

.994198 

.995579 

.997046 

.041859 

.989455 

.992289 

.994109 

.996107 

.055714 

.985974 

.989744 

.992171 

.994793 

.074155 

.981379 

.986356 

.989582 

.993026 

.098701 

.975260 

.981871 

.986122 

.990741 

.131371 

.967202 

• 975905 

.981549 

.987666 

• 174854 

.956595 

.967981 

.974281 

.983536 

.232731 

.942654 

.957496 

.967353 

.978117 

.309765 

.924480 

.943676 

.956601 

.970900 

.412298 

.900963 

.925521 

.977487 

.961324 

.548769 

.870840 

.901780 

.923601 

.948489 

.730410 

.832991 

.871117 

.898954 

.931629 

.972176 

.786407 

.832080 

•866865 

.909231 

1.293971 

.730797 

.783488 

.825723 

.879799 

1.722267 

•666936 

.724985 

.774219 

.841413 

2.292342 

.597051 

.657693 

.712129 

.792279 

3.051106 

.524404 

.584387 

.641035 

.731374 

4.061021 

.452652 

.509051 

.564575 

.659777 

5.405230 

.384325 

.435216 

•486945 

.581021 

7. 194348 

.320979 

.365417 

.411695 

.499720 

9.575701 

.263678 

.301455 

.341523 

.420331 

12.745267 

.213191 

.244575 

.278317 

.346376 

16.963960 

.169876 

.195418 

.223168 

.280303 

22.578970 

.133684 

.154127 

.176499 

.223239 

30.052592 

.104149 

.124756 

. 138040 

.175498 

40.000000 

.080462 

.093050 

.106967 

.136563 
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TABLE I CONTINUED 


TRANSMISSION PROBABILITIES FOR CYLINDRICAL DUCTS 
AT VARIOUS SPEED RATIOS 

(CALCULATED FROM UNIVERSITY OF TORONTO REPORT NO. 88» 
DECEMBER 1962* BY J. H. DELEEUW AND D. E. ROTHE) 


L/A 


K ( S » L / A ) 



S=4 • 0 

S = 5 • 0 

S=10 

S-20 

.010021 

.999324 

.999422 

.999728 

.999854 

.013338 

.999078 

.999232 

.999621 

.999817 

.017752 

.998780 

.999003 

.999504 

.999755 

.023628 

.998319 

.998668 

.999334 

.999668 

.031449 

.997809 

.998212 

.999126 

.999556 

.041859 

.997074 

.997655 

.998812 

.999407 

.055714 

.996051 

.996861 

.998422 

.999218 

.074155 

.994749 

.995837 

.997901 

.998948 

.098701 

.993071 

.994415 

.997208 

.998606 

.131371 

.990710 

.992590 

.996293 

.998146 

• 174854 

.987687 

.990126 

.995078 

.997530 

.232731 

.983614 

.986862 

.993430 

.996716 

.309765 

.978170 

.982532 

.991253 

.995628 

.412298 

.970961 

.976744 

.988376 

.994182 

.548769 

.961334 

.969042 

.984521 

.992262 

.730410 

.948588 

.958846 

.979401 

.989701 

.972176 

.931699 

.945256 

.972597 

.986288 

1.293971 

.909277 

.927241 

.963531 

.981747 

1.722267 

.879753 

.903438 

.951487 

.978051 

2.292342 

.841134 

.872023 

.935498 

.967719 

3.051106 

.791366 

.830957 

.914288 

.956995 

4.061021 

• 729290 

.778163 

.886246 

.942842 

5.405230 

.655847 

.712652 

.849368 

.923969 

7.194348 

.574787 

.635915 

.801390 

.899044 

9.575701 

.491180 

.552293 

.740265 

.866123 

12.745267 

.410023 

.467380 

.665846 

.822983 

16.963960 

.335220 

.386354 

.581363 

.767235 

22.578970 

.269265 

.313020 

.493029 

.697465 

30.052592 

.213006 

.249408 

.407378 

.615228 

40.000000 

.166509 

.195940 

.329503 

.526056 
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C TRANSMISSION PROR AB I L I T I F S FOR CYLINDRICAL DUCT c 

COMMON I S TA RT ( ? ) 

OTMFNSION COI !NT ( 8 , 7 ^ ) 

1 FORMAT ( F0.4,5XF9.4»SXF0.4,SXOF.O) 

2 FORMAT ( FO.4 ,FXF9.4,5XE0.4 ,5XF9.4 ) 

3 FORMAT ( 1X°H COLLISION FX9HFXIT AT L 7X Q HFX I T AT o7Y]nHFxiT AT L?) 

4 FORMAT (2X13, 12XF5.0, 12XFO.P, 12XF5.0) 

5 FORMAT ( nX3HX=L?XF5.O,8X3HX=02XF9.0,8X4HX=L22XF5.P*7X3H0UT2XFS.p) 

6 FORMAT (014.7) 

7 FORMAT ( ?H?7,?XlHT?yoi 4 . 7 . 5X 1 H'/2X r 1 4 • 7 , S X 1 H'-'7X C 1 4 .7 ) 

8 FORMAT r?H70,?XlHT2X c 14.7,FXlMV?xri4,7,^VlHW?Xoi4.7) 

o Fop M A T ()u ) 

77 FORMAT (FMX(IC), 7X014.7, 5XRHC0LL NU'S 2X18) 

151 FORMAT (05.3) 

152 FORMAT ( 3HL1=F9.4,5X3HL2=F9.4 ,5X2HS=F9.4, 5X5HTOT AL7XF5 . 0 ) 

153 FORMAT { 1x7HA = FO«4,6X , HR = F9.4,Fx?0 | C = F o .4,5X?H0 = F o ,4 ) 

260 FORMAT ( 1 X 4 4 H M 0 N T F CARLO SOLUTIONS OF FRFF M 0 L F C U L A P FLOW) 

761 FORMAT (1X2 5HTHROUOH CYLINDRICAL OIICTS) 

76 7 cnoMAi ( 1 yl 6 mamo cdf^p ratTPC) 

760 FF'P-'AT ( lyl r-HPO CDFFp R^TIOC) 

770 FOiDV at ( 1 U1 ) 

771 FORMAT (14) 

27? FORMAT ( 1 X7 5HSPFF0 RATIO = K ( 1 . 1 *4 700 ) ) 

604 FORMAT ( 6 H X ( I C 1 ) * F 1 /, . 7 , 6HX ( I C 7 ) , F i 4 . 7 , 8 HCOL L MUM, 13) 

6 05 FORMAT ( ! WYFI 4 .7 , 1 H70i/ f4 7 , 1 HXF1 4. 7 »?H01 oi 4.7 , ?H0?P1 4 . 7) 

606 Format ( 7MU1 F) 4 .-7 , 7 H 1 1 7 F 1 4.7,7Hi!8ri4.7,7MAAPlA.7,7HAO0]4.7,Fl4 # 7) 
61] FORMAT (FHX(IO), 2X014.7, 7HC0LL NO, 7X1°) 

6^0 FORMAT (OP, P X o A ) 

P700 FORMAT ( F ° • 4 ) 

R 7 o] FORMAT ( 1 X76HAMOLF 00 ATTACK T N OFopr^e;, cq,4) 

PI = 3. 141 597 A 536 
150 XL 2 T = 0. 

7800 FORMAT ( 5X6HI START , 5X04, 5X08) 

XOT = 0, 
x L T = 0. 

OUT = o, 

A T O T = O, 

ROAD 1 , O, OL 2 , SI , T M A X 

RFAP 7, A, B , CL7, DL7 

RFAD 650, ISTART(l), ISTART(7) 

ROAD 771, MCOL.L 
RFAD 870A , ANOLF 

O P I M T 7 7 0 
PR I M T 7 AO 
D P I M T 7 A 1 

IF (SI) 765,766,265 


FIGURF 7 


FORTRAN PROGRAM FOR momjf CARLO CALCULATIONS OF 
FRFE MOLECULAR FLOW THROUGH CYLINDRICAL OUCTS 
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26 ^ 


266 
2 69 


500 

2010 

20 ]]. 

201 ? 

2 n l * 
21 4 p 

2555 


pssp 


7 5 50 
2 5 57 


2141 
2 5 56 


poop 


PPIMT 26 7 
PPIMT 27? 

GO TO ?so 
D P I M T ? A S 
PPIMT o 
PPIMT 9 

PRINT 152, Q, QL 2 * SI, TVLAX 
PRINT 159, A, R, CL? , PL? 

PRINT 8701 , A N G L F 
d p I w T 9 

PO 500 J = 1,9 
no 500 < = 1 , NCPLL 
COUNT (JfK) = 0 # 

A A T A C = A MG L F * • 0 1 7 4 5 3 2 0 2 5 
VI = { 1 • - 2.*FLAT ( ISTART ) )**A 

71 = (1. - ?.*FLAT( ISTART) >*A 

0 = Y 1 *Y 1 + 71*71 

1 COLL = 1 

if (o - a*a ) ?oii ,?om,poio 
R 1 = FLATtI START) 

R? = FLAT (ISTART) 

I F ( R 1 - R?) 2013,2.01 1 ,201? 

CT H r T = PI 
GO. TO 2142 
CTHFT = 07 

S T H F T = SQRTFM. - C TH r T **? ) 

sn = 1 # 

IF (SI) ^55,2556*255^ 

XT = 0. 

YT = 0. 

7 T = o, 

002557 J = 1*1? 

XA = FL AT ( I S T ARJ ) 

YA =. FLAT ( j START ) 

YT = YT + va 
7 A = FLAT ( I P T A P T ) 

7 T = 7 T + 7 A 

XT = XT + XA 

V X = XT - 6* 

VY = vt ~ 6. 

V 7 = 7T — 6* 

SP = SORT ( V X * V X + \/Y*'/Y + \'7*\/7) 

I r ( F L A T ( T S T A P T ) - .5) 2 1 4 1 * 2 0 1 1 , 2 * r > 5 

CTHFT = “TTHPT 

PS I = F L A T ( ISTART)*2.*PI 

AA = CTMrjrr FP + S I * 1 . 1 54 7 OQ*C OS F ( A A T A C ) 

AP = STHFT*CQSF (PSI )*SP + SI*1 . 1 5470O*s i NF ( a.AtAC) 

AC = STHFT * SlNFfPs r )*SP 

OF NO = SORT(AA*AA + AP*.AP + AC#4C) 

1.11 = A A /nrM o 

IF nn - n. ) ? 0 1 ] ,201 1 ,ooop 

•JP r A^/pPMO 
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FTPURF ^ CONTINUE 



U3 = AC/DFMO 
AT. ST = Atct + 1. 

20? 1 HFLTA = H2/M1 
DFLTB = MVtll 
?n ?0 XI = n. 
x i x l = c 
xixo = o. 

V I X L = DFLTA*0 + Y1 
ZIXL = DFi.TR *0 + 71 

02A IF ( Y I X L * Y I XL + 7IXL*7IXL - 1.) A 1 no , A 1 00 * 302 6 
3026 T = DF L T A * Y 1 + DFLTR*71 

V = DFLTA*D C LTA + DFL T R #nFLTB 
W = Y1 *Y1 + 71*71 - 1 . 

3029 P = SQRTF ( T*T - V#H) 

XIC = (-T + P ) /V + -XI 
3 0 3 A IF (XIC - Q) 3035,3035 *3036 
3036 PRINT 37, XIC , ICOLL 
GO TO 7115 

300 3 YIC = DF[.t A* ( X IC - XI) + Y 1 
Z I C = DFLTB* (XIC - XI ) + 71 
3039 STAU = 7 1 C 
C T Alt = Yjr 

3025 ICOLL = ICOLL + 1 

IF (ICOLL - MCOLL) 3 3 A 1 ,7115*71 l c 
3A1 IF (XIC ) 3036,3036, OOAO 

OnAO R l = F! AT ( I START ) 

R2 = FIAT (I START) 

IF ( R 1 - R 7 ) 3042 , 3RA0 , 30A1. 

3 0 A 1 CTHFT = R 1 
GO TO 30A 3 
30A2 CTHFT = R2 

3 0 A3 STHFT = S0RTF(1. - CTHFT**2) 

3 0 A A PSI = FLAT( ISTART)*?»* D I 
IJ1 = S I M F (PSI ) * S T H F T 

I I? = -CTHFT*CTAU - S T H f T * C 0 S F ( R S I ) * S T A ( i 
U3 = -CTHFT *ST A IJ + STHF T*COSF (PSI ) *CTAIJ 
DFLTA = H2/U1 
DFLTB = i '3/1)1 
IF (1)1 ) 3060, 30A5, 3045 

30A5 YIXL ='DFLTA*(Q - XIC) + YIC 
ZIXL = D r L T B * ( 0 - XIC) + 7 T C 
Y1 = YIC 
Y 1 = YIC 

71 = 7 I C 

1052 IF ( YI XL*YI XL + 7TXL*7TXL - 1.) A 1 0 0 , A 1 n o , 3 06 P 
3060 YIXO = DFLT A* ( -X IC ) + YIC 
ZIXO = 0 f (_ t R * ( - X I C ) + 7IC 

XIXO = 0. 

XI = XIC 
Yl = YIC 
71 = Z IF 

1063 I F ( YI X0* v I XO + 71X0*71X0 - 1.) 5 1 1 0 , S 1) n , SO 6 p 


CIGURF 3 CONTI Nl.iFO 


?068 XIXO = 0. 

3069 T = 0FLTA*Y1 + PFLTR*71 

v = DFLT a*of| t a + PFLT Q * r ' p LT^ 5 
IF (IJ1) 8691 ,*039, 860? 

3691 XIC = -APFF( <?.*T)/V) + XI 
GO TO 3034 

360? XIC = A R 9 F ( ( ? • * T ) / V ) + XI 
GO TO 30 o 4 

100 IF (YIXL*YIXL + 7IXL*7IXL - B*R) 4101,4799,4799 
410] COUNT ( 1 ♦ I COLL 1 = COUNT ( 1 , 1 COL L) + 1. 

X| T = XL T + 1, 

IF (OL?) 71?o, 71?o, 6170 

4709 I COLL = TOOL 1 + 1 

IF ( ICOLL - UCOLL) 48 00 » 7 1 1 5 , 7 11 c 

4800 R 1 = FLAT(ISTART) 

R? = FLAT(IFTAPT) 

I F ( R 1 - P? ) 4801 ,4800,480? 

4801 CTHRT = - p? 

GO TO 4800 

4807 CTHFT = - PI 

4808 STHFT = c ortf ( 1 . - CTH C T**? ) 

P.8 1 = FLAT! I9TAPTJ*2.#DI 

111 = CTHFT 

1.1? = 8THFT#C0FF { PF I ) 

U 3 = 9 T H f T * 8 I M F ( P 8 I ) 

OFL TA = M 7 / I 1 1 
OFLTP = " 8 /i»l 

4804 Y I XO = OOLTA^j.xjyL) + Y I XL 
7 1X0 = P r l. T9* (-X I XL ) + 7 1 XI 
XI = XI XL 

Y 1 = Y I XL. 

7 1 = 7 I X|_ 

IF ( Y I X 0 « y I X 0 + .7 1X0*7 1X0 - 1.) 8 110,81 1 0,4808 

40 OF T = OFl.T A*Y I XL + OF L. TO* .7 I XL 

V = OF IT A *OFL T A + of| T°*oFLTP 
= YIXL* V IXL + 7 I x L * 7 I y L - 1, 

XIC = ( -T — 90PTF ( T*T - V*H))/V + XI. XL 
GO TO 3086 

>110 IF (YIXO*YlXO + 7 1X0*71X0 - A * A ) 8 1 1 1 , 8 p 1 0 , p p } n 
5111 COUNT (2, 1 COLL) = COUNT ( ? , T COLL ) + 1. 

X OT = XOT + 1 . 

0,0 TO 718 0 

6810 I fo l t = TOOL | + 1 

IF (TCOLL. - MCOLL) 5 8 11,7118,7116 
^811 Rl = FL AT ( I ot.APT ) 

R 2 = FL AT ( I 8 T A R T ) 

IF (Ri - op ) 6 8 1 ? , c 8 1 1 , 8 P 1 8 

K 817 CTHFT = P7 

Go to 8 bi 4 
8818 CTHFT = PI 

coi A 8 T H r T = cnojF ( i, _ rruFT,#?) 
p 8 1 - FI A T ( I <T APT ) *? . *P I 


F ICL'PF 


7 


CONTI N 1. 1 F 0 



ill = r t m r t 

1 1 7 - c t i-’ c T c o . r r ( r> t ) 

i n = c j u r t c ] * 1 cr ( D c j ) 

Y 1 = YTY^ 

7 ] = 7 J Y O 

C,r\ j r\ ? A 7 1 

617^ v L 7 = Y I y f_ 4 0( 7 

Y Y | 9 = or I TMMXL? - YIYI ) 4 v T v I. 

7 V ( 7 = "> r I T Q * ( Y I 7 - Y J V \ ) + 7 T Y I 

T r ( v y | ? -v. v y j p + 7 v| 0 7 v | 9 _ r\ 0 -y. r\ "> ) «QCA^m 

6 171 jrA| | = TrA) | + 1 

I r (Trni_| - M r 0 | | ) M 7 9 , 7 1 1 r, > 7 1 1 r. 

6179 T = ^FL T a ->Y T YL + P r LT‘^7 I YL 

V = P c L T A * n F l T A + n r i T"-nr| T " 

= Y I XL * Y I X L + 7 T XL *7 I x L - n ?*n 7 

O = CA 7 T^{T^ T - \t 4 11 ) 

X I c 9 = ( ( -T+P ) /VM + v T y | 

y i — v t yi 

YI r YTYI 
7 1 = 7 T Y I 

6 7?^ Yir? = n r l TA^(Yir? - Y 1 ) 4 YI 

7 Tr? = nr| T p^ (V jr? - yi ) 4 71 
c j a 1 1 = 7 r r ? / r i ? 
c t a !. ■ = Yir. 7 /rL? 

617 9 PI = F L A T ( I c T A P T ) 

09 = o | A T ( | C T a o J | 

T r (Pi - 09) 6 i 7 c * * 1 7 ^ * '■ 1 ^ /i 

6174 CTH r T = oi 
in 6176 
6 1 7 R rjH-T = P9 

61 7 £ PJH^T = r n n T p ( ] # - r T tiF T **? ) 

d c I = Fl a T ( T’C j apt ) * 7 . #P I 

in = c J M r ( p r - J ) # r T M c T 

II 7 - «. r J • ! r T T A I ' - c yu r j 4 a n 7 r ( n c. j 1 4 c T M 1 

117 = - c T I • 1 r T K C T a 1 I 4 r T U r T 4- A A r <T ( n C t 1 .>? A T A 1 I 

O r | -J- a = ' i 9 / I 1 1 

or | rn - no/m 
IF ( i *1 ) f 1 7 7 , 6 1 7 9 , 6 1 7 9 
6 17 7 YYF ,IM = nr ! I A* (Yiy[ - v ] r ? } 4 Y I T 9 

7 v o m = nrj T n * ( Y T y L. - ytG 2) + 7ir^ 

IF ( v y a m # yy o *■ 4 7 y 0 v 4 7 y a v - rt ?*r\ 7 ) 6 1 40»6i op, 6P1 l 

6179 YY|° = P - IT A 4 ( Y |9 - Y J r 9 ) 4 Y T T 9 

7 Y L 9 = P F t. T ° 4(Yf9 _ Y T 9 9 ) 4 7 T 6 P 

I r ( V y [7 4 V Y |_ 7 4 7 y L_ ? -«- ? v | 9 - C \ 7 * r \ 9 ) 6 0 r > 0 , a p.i 1 , 6 0 *» 1 

■ IPO [F ( YXCM^vynM 4 7 X M 7 Y -P A’ “ p n ) ^1^1 ,61^1 *6 1^9 

6 1 Q 1 Y T X l = Y X 0 N* 

7 T Y 1 = 7Y9M 

m jo a poh 

6 IR? I r OL I = T r ° L I + 1 

jr (jrrj i - trn [ I ) 6 1 9 9 ,711 K ,71 1 r 

61P9 Pi - f| at ( I c T A P T ) 

P? = CT| at ( I CJAPT ) 


ftgiup- 9 rowj t m» i r n 



IF ( R.l - R? ) 61. 84,61 P”’ ,61 R8 

f,m f T H F J = ”1 

00 r A A 1 P 6 
AIR 4 CTH^T = » ? 

• 6 1 8 A 8TH P T = F0RTF(1. - CTH p T**?) 

PSI = C L^T( I C .TART)*?.*PI 

t.n = cth c t 

! I ? = FTwr T *coRF ( P p I ) 

!J? = STHFT*SI N| F ( P p I) 
nr|_TA = n?/in 

D r | T n = " R / 1 n 

YIXL = YVO'i 
7 1X1. = 7 Y o N 
AO TO 6170 
ART] XI = Y T A 7 
Y 1 = YIC7 

7 1 = 7 I C? 

1 ro| L = T o o | | + 1 

I p (iron - moou ) 6810,71 i R,7n R 

6810 T = DFL. T A*Yl + D e I.TP*71 

V = 0 F L T A *D F L T A + D r LT p *PFLTP 

IF ( HI ) 68.71 ,7.1 15, AR7? 

68.? 1 XIC? = -A 0 8F( (?.*T )/\f) + XI 
AO JO 6770 

68 ? 7 XIC? = An e.F ( ( ?.#T ) /V) + XI 
AO to 6770 

A 8 c o if ( Y X L ?* v v L ? + 7XL?*7XL? - PL7#n|?) A8A 0 , A 8 8 1 , A 8 R 1 
A R 5 1 IOOLL = I r o L L + 1 

IF (IOOLL - MCOLL) 688?, 71 15,711 8 
6887 R] = FLAT ( I5TAPT) 

P? = FLAT ( I START) 

IF (Pi _ d 7 ) A887,68 c .7,AR r -A 

68 8 R ATH r T = °7 
oo to a^RR 
F R R /. f T H c T = oi 

A 8 8 8 c T H r T = c OR T F (1 . - CTWfT**-?) 
pci = FLAT( ICTA.RT)*2.*PI 
'll = - C T h F T 
M? = CTH p T*CO c F ( pc I ) 

HR = 8 T FI 0 T * 8 i m F ( P S I ) 

OFLTA = "?/ll] 
o r l T P = i ' 8 / 1 ! 1 

Y YON = O r | j A * ( Y I X L - Y I 7 ) + YXI ? 

7 YON = n r |T°'*(XIXL - XI 2) + 7 XL. 7 

I F ( YXO w *Y YON + 7 X C M * 7 Y O N - CL 2 *01 ? ) 6 1 8 0 , 6 1 8 0 * A 8 c A 

A R 5 A T = OFLTA #v XL? +0FLT p *7XL? 

V = noLTA^OFLTA + PFLT n #P r LT n 

'■! = YXL ?*YX L? + 7XL?*7XL? - Cl.7*CI. 7 
XIC? = (-T-SORTF (T*#? - V* '•' ) ) / V + XL? 

71 = 7 x|, 7 

Y 1 = Y x I 7 

XI = XI? 
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FI-A'iRF 8 rONTTMt.'Fn 



60 

7115 
7170 
77 1] 

7140 

9np,fi 


GO TO 67 n o 

COUNT (7, t COL I) = roimj (fi^iroj | ) + i, 

VI 7T = XL7T -f 1. 

GO TO 7170 

OUT = O' IT + 1 • 

IF ( A T c T - TV AX) 7^1 o,77i 1 *7 ? ^ 1 
p p I m t 7 

00714-0 J = 1 9 MfOIJ 

< = J - 1. 

PRINT 4, fOtiMT(T,J), roiMT(^j) 

P D | M T 0 

op I N'T K , Y[.T, Y 0 T , Y| i OMT 

po I M T o 

POINT 7^n, I c TAPT(l). TCTAOT(O) 

GO TO l ^ n 
r m p 


F t Gi I or ^ r o m t I Ku * rn 
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MONTE CARLO SOLUTIONS OF FREE MOLfcCULAR FLOW 

THROUGH CYLINDRICAL DUCTS 

AND SPEED RATIOS 

SPEED RATIO = Ml. 1547Q0) 


1= 2.0000 

12= 0 S= 

.4640 

TOTAL 

9999 

A« 1.0000 

8 = 1.0000 C = 

0 

D = 

0 

ANGLE OF ATTACK 

IN DECREES 40.00 0 0 





COLLISION 

EXIT AT L 

EXIT AT O 

EXIT AT L2 

0 

2203 

0 

0 

1 

1209 

2126 

0 

2 

853 

975 

0 

3 

533 

542 

0 

4 

319 

320 

0 

5 

178 

180 

0 

6 

98 

118 

0 

7 

70 

74 

0 

6 

42 

45 

0 

9 

25 

17 

0 

10 

15 

11 

0 

11 

11 

6 

0 

12 

7 

9 

0 

13 

2 

3 

0 

14 

1 

3 

0 

15 

1 

1 

0 

16 

1 

0 

0 

17 

0 

0 

0 

18 

0 

0 

0 

19 

0 

0 

0 

20 

1 

0 

0 

21 

0 

0 

0 

22 

0 

0 

0 

23 

0 

0 

0 

24 

0 

0 

0 

25 

0 

0 

0 

26 

0 

0 

0 

27 

0 

0 

0 

28 

0 

0 

0 

29 

0 

0 

0 

30 

0 

0 

0 

31 

0 

0 

0 

32 

0 

0 

0 

33 

0 

0 

0 

34 

0 

0 

0 

35 

0 

0 

0 

36 

0 

0 

0 

37 

0 

0 

0 

38 

0 

0 

0 

39 

0 

0 

0 

40 

0 

0 

0 

41 

0 

0 

0 

42 

0 

0 

0 

43 

0 

0 

0 

44 

u 

0 

0 

45 

0 

0 

0 

46 

0 

0 

0 

47 

0 

0 

0 

48 

0 

0 

0 

49 

0 

X=L 5569 

0 

X*Q 4430 

0 

X = L2 


OUT 
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FIGURE 4 - TYPICAL COMPUTER RESULT 



Number of Molecules XI0~ 
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FIG. 5. DISTRIBUTION OF SPEEDS FROM MONTE CARLO CALCULATION 







Theory 
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FIGURE 7. TRANSMISSION PROBABILITES FOR 
CYLINDRICAL DUCTS AS A FUNCTION OF SPEED RATIO 



Calculated from UTIA Report # 88 
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L/A 

FIG. 8. TRANSMISSION PROBABILITIES FOR CYLINDRICAL DUCTS AS fi 
FUNCTION OF LENGTH TO RADIUS RATIO, L/A, WITH SPEED RATIOS 



Calculated 
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FIG. 9. TRANSMISSION PROBABILITIES FOR CYLINDRICAL DUCTS AS A 
FUNCTION OF LENGTH TO RADIUS RATIO, L/A, WITH SPEED RATIOS 
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FIGURE 10. TRANSMISSION PROB ABILITES FOR 
CYLINDRICAL DUCTS AS A FUNCTION OF SPEED RATIO 



UTIA 

Monte Carlo 
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FIGURE II. TRANSMISSION PROBABILITES FOR 
CYLINDRICAL DUCTS AS A FUNCTION OF SPEED RATIO 



UTIA 

Monte Corlo 
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FIGURE 12. TRANSMISSION PROBABILITES FOR 
CYLINDRICAL DUCTS AS A FUNCTION OF SPEED RATIO 



UTIA 

Monte Carlo 
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FIGURE 13. TRANSMISSION PROBABILITES FOR 
CYLINDRICAL DUCTS AS A FUNCTION OF SPEED RATIO 




UTIA 

Monte Carlo 



FIGURE 14. TRANSMISSION PROBABILITES FOR 
CYLINDRICAL DUCTS AS A FUNCTION OF SPEED RATIO 





UTIA 

Monte Carlo 
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FIGURE 15. TRANSMISSION PROBABILITES FOR 
CYLINDRICAL DUCTS AS A FUNCTION OF SPEED RATIO 





FIG. 16 . TRANSMISSION PROBABILITIES FOR CYLINDRICAL DUCT 
AT A CONSTANT SPEED RATIO FOR VARIOUS ANGLES OF ATTACK 
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FIG. 17. TRANSMISSION PROBABILITIES FOR ORIFICE RESTRICTED CYLINDRICAL TUBE 
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FIG. 18. TRANSMISSION PROBABILITIES FOR TWO CYLINDRICAL DUCTS IN SERIES 




FIG. 19. TRANSMISSION PROBABILITIES FOR TWO DUCTS IN SERIES AT VARIOUS SPEED RATIOS 
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